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Abstract 

Dynamical Ensemble Equivalence between hydro dynamic dissipa- 
tive equations and suitable time-reversible dynamical systems has 
been investigated in a class of dynamical systems for turbulence. 
The reversible dynamics is obtained from the original dissipative equa- 
tions by imposing a global constraint. 

We find that, by increasing the input energy, the system changes from 
an equilibrium state to a non-equilibrium stationary state in which 
an energy cascade, with the same statistical properties of the original 
system, is clearly detected. 

PACS number 47.27.Jv, 47.90. +a, 05.45. +b 

To the memory of Giovanni Paladin. 
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1 Introduction 



One of the most important open problem in classical physics is the under- 
standing the statistical features of fully developed turbulence (FDT). 

A fully developed turbulent flow is a dissipative system described by the 
Navier-Stokes (NS) equations in the limit of high Reynolds numbers (Re). 

On one hand, direct numerical simulations of turbulent flows are strongly 
limited due to the huge amount of excited degrees of freedom: a simple 
argument due to Landau, shows that the number of degrees of freedom, 
which should be taken into account for a correct description of a turbulent 
flow, increases as Re 9 / 4 . 

On the other hand, analytical attempts to derive the multi-points ve- 
locity probability distribution function have repeatedly failed due to the 
strong coupling regime and due to the highly non-gaussian probability dis- 
tribution functions (PDF) developed at small scales by the velocity field [p]]. 
Phenomenological approaches, or simplified dynamical and deterministic sys- 
tems, have been therefore often used for studying the mechanisms generating 
the turbulent energy cascade. 

From the analytical point of view, the main obstacle to the possibility of 
performing the 0-th step toward a theory of turbulence is certainly connected 
to the strong dissipative and far-from-equilibrium character of 3d turbulent 
flows. In 2d turbulence, where energy is almost not dissipated at all, some 
analytical tools based on quasi-equilibrium statistical ansatz have indeed 
been developed [[|. 

Strongly chaotic dynamical systems as the Anosov systems are the only 
cases where, although being still dissipative and chaotic, analytical tools have 
been developed with relative success, at least in the case of low-dimensions 

I- 

Recently ||, Gallavotti and Cohen proposed that a chaotic high- dimensional 
dynamical system in a stationary state can be regarded as a smooth dy- 
namical system with a transitive Axiom-A global attractor or, if it is time- 
reversible, as a smooth transitive Anosov system, as far as macroscopic prop- 
erties are concerned. This is the so-called Chaotic Hypothesis. 

This hypothesis, as the ergodic hypothesis, can be proved only in very 
particular systems. Nevertheless, it is interesting to analyze some of its 
consequences. 

Likewise, Gallavotti [§, U conjectured a Dynamical Ensemble Equiva- 
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lence between some dissipative systems (in this paper we will refer only to 
3d Navier-Stokes equations) and their non- equilibrium but time-reversible for- 
mulation. Here, equivalent must be meant that the averages of local variables 
in the two systems, i.e. the original one and its time-reversible formulation, 
are the same, in a suitable limit. For hydrodynamic dissipative systems this 
limit is that of FDT, i.e. Re — > oo. 

The reversible dynamics is obtained from the original dissipative equa- 
tions by imposing a constraint, such as to keep constant in time those macro- 
scopic observables (as the total energy) which would have only had stationary 
averages in the original systems. 

Having a time-reversible system and applying the chaotic hypothesis, 
some large- deviation properties of the fluctuation of the entropy-production 
rate in the system can be proved [|J. 

In Navier-Stokes equations, the reversible dynamics is achieved by in- 
troducing a sort of "eddy viscosity" which removes the input energy with 
perfect efficiency Viscosity becomes non-positive defined and strongly cor- 
related with the large-scale flow where energy is injected. 

As far as we know, the idea of reversible NS equations was introduced 
for the first time by She and Jackson 0. They did not exploit global con- 
straints, but imposed that the energy contained in each "momentum shell" 
was constant. 

Let us also mention that a constraint of constant energy has been im- 
plemented by using the Kraichnan's eddy-viscosity parameterization || in 
low-resolution large-eddies simulations of NS equations in ||. With such a 
parameterization of the viscosity one has very weak fluctuations of the energy 
(less than lof high-resolution numerical simulations. 

From our point of view, the interest in models with global constraints 
stems from the possibility to describe a global macroscopic dissipative and 
irreversible physics starting from a deterministic reversible dynamics. The 
approach can be seen as a bridge from microscopic reversible dynamics to 
macroscopic irreversible dynamics and, more interesting, a possible system- 
atic tool for going with continuity from a pure-equilibrium and conservative 
systems to a strong dissipative and far-from-equilibrium time evolution. 

In this paper, we investigate these ideas in Shell Models, i.e. a class of 
simplified dynamical systems for turbulence (for a recent review see JTUJ, for 



a tutorial introduction see [fLl]| ). In particular, we will analyse in details 



the smooth transition from the equilibrium system at zero viscosity and zero 
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external forcing to a (formally)-reversible systems which possess anyway a 
non-equilibrium flux of energy from large to small scales. 

The paper is organised as follows. In section 2 we briefly review the sta- 
tistical mechanics of a perfect fluid and the ideas presented by Gallavotti in 
(!) H § concentrating only to the case of Hydrodynamical systems (Navier- 
Stokes eqs.). In section 3 we discuss Shell Models philosophy and we describe 
Gallavotti's implementation to our case. In section 4 we present our numer- 
ical results. Conclusions follow in section 5. 

2 Equilibrium and Non-equilibrium Statisti- 
cal Mechanics 

In a 3D perfect fluid, i.e. with vanishing external forcing v = 0, the evolu- 
tion of the velocity field is given by the Euler equations which conserve two 
quadratic functional, the kinetic energy and the helicity: 

E= l -<v 2 > H=~<vu>, (1) 

where uj = V x v is the vorticity. In this case, it is possible to construct 
a statistical mechanics as for a gas: by using the conservation laws and the 
conservation of the volume in phase space one obtains a gaussian distribution. 

For simplicity let us start by neglecting the helicity conservation. To 
be explicit, let us consider an incompressible inviscid fluid in a cube with 
periodic boundary conditions, so that the velocity field can be expanded in 
Fourier series as 

^■(x) = L- 3 / 2 ^e ikx ^(k) (2) 

k 

with k = 2nn/L and n = (m, n%, TI3), where n, are integers. The variables 
Vj(k) are not completely independent, since from the incompressibility con- 
dition and the fact that v(x) is real, it follows that 

k-v(k) = and v(k) = v*(-k) 

In any case, it is straightforward to introduce a new set of independent 
variables X a , where now a labels the spatial component and the wave vector. 
By using an ultraviolet truncation, v(k) = for k > k mSuX , and by introducing 
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(0) in the Euler equations one obtains a set of ordinary differential equations 
(ODEs) with the structure 

— TT = M a bc X b X c 

dt Tl. 

b,c 

where M abc = M acb and M abc + M bca + M cab = with a = 1, • ■ • , N ~ k z max . 
We stress the fact that the ultraviolet truncation is necessary in order to 
avoid the infinite energy problems of classical field theory. 

It is easy to verify that [| preserves the volume in the phase space as well 
as the energy, namely 

£^(^1=0 and ^ = -^E^ = 
V dX a \ dt j dt 2dt^ a 

These conservation laws are sufficient to construct the probability distribu- 
tion of the variables {X a } J0|: using the ergodic hypothesis, one obtains the 
microcanonical probability measure 



It is well known that, in the limit iV — > oo, this is equivalent to the canonical 
measure 

P c ({X a }) ~exp 



where the Lagrange multiplier j3 satisfies the relation 

In two dimensions, the helicity H = and there exists a second conserved 
quantity, the enstrophy 

n = l fcu 2 d 2 x 



2 

which is the mean square vorticity. In terms of the X variables, it can be 
written as 

n = \Y,klxl 
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As a consequence, the microcanonical probability measure in 2d is 

and the corresponding canonical measure is 

P c ({X a }) ~ exp - (| E ^ + y E ^) 
where the Lagrange multipliers satisfy the relation 

< Xl >-- 



In 3-d one can repeat a similar argument, taking into account the helicity 
conservation. In this case, being the helicity non positive defined, one as to 
assume suitable constraints for the generalized temperature related to the 
helicity @. 



The above results are, both in 2d and in 3d, well reproduced by numerical 



simulations [12]. 



The limit v — > (equivalent to Re — > 00) is singular and cannot be 
interchanged with the limit N — > 00. Therefore, the statistical mechanics of 
an inviscid fluid has a quite limited relevance on the behaviour of the Navier- 
Stokes equations at high Reynolds number. Recently some authors proposed 
the use of conservative statistical mechanics to justify some behaviours of real 
fluids, e.g the Jupiter's red spot and the emergence of organized structures 



14], |T|, [Lq| . The applicability of this approach is limited to some particular 
quasi equilibrium two-dimensional situations. 

On the other hand, both from phenomenological arguments and exper- 
imental results, we know that the statistical mechanics of fully developed 
turbulence has peculiarities rather different from the usual statistical me- 
chanics of conservative systems. In the limit of FDT the energy fluctuates 
around its mean value and in addition one has an energy cascade from large 
to small scales. 

The turbulence is described by a dissipative system (essentially a high 
dimensional truncation of the Navier-Stokes equations with |k| < k max = 
0(i?e 3//4 ) in which the volume in the phase space is not conserved. Let us 
stress again that the two limits v — > and k max — > 00 must be take in a 
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suitable way in order to obtain the correct physical result for the turbulence. 
If one wants that the mean energy dissipation is 0(1) in the limit Re — > oo 
one has to take k max > 0(Re 3 ^ 4 ). 

In order to have a statistical stationary state one needs two basic ingre- 
dients: a 'friction' mechanism and a coupling with an external forcing or 
'reservoire'. A typical example of statistical stationary state is given by con- 
ductive systems where an external electric field and a friction mechanism, 
mimicking the electrical resistivity, leads to a macroscopic steady current. 

Recently time-reversible and conservative systems have been introduced 
in the issue of the non-equilibrium statistical mechanics of stationary state. 
For sake of self consistency we recall one of the simplest system of this class. 
Let us consider N independent particles of mass m, with coordinates and 
momenta and respectively, on a square domain (in 2 or 3 dimensions) 
with periodic boundary condition for the variables q^. Using a suitable ex- 
ternal potential V^(q) we can mimic the elastic scattering with rigid obstacles 
in order to have basically a "Lorentz gas". The equation of motion are: 

dqi 1 
dpi = _dV_ 

dt dq { 1 ] 

The system is chaotic and one can expect the usual microcanonical distribu- 
tion. Of course there is no net current. In order to have a current in the x 
direction it is necessary to add in the eq. (f|) a term Eei, where ei is the 
versor in the x direction. At the same time, if one wants to focus on sta- 
tionary aspects, some energy-loosing mechanism must be added. Standard 
phenomenology would suggest the insertion of a viscous irreversible term of 
the form — api in the equation of motion governing the evolution of momenta. 
In this way, one is naturally lead to a stationary state. 

Recently, in [ 171 , the idea of mimicking this behaviour by means of an ex- 



actly conservative and reversible physics has been proposed by using instead 
of a constant viscous coefficient a, a perfect energy-sink, correlated with all 
scales and able to reabsorb instantaneously all excess of energy injected by 
the forcing term in the system. This ideal viscosity must acquire an explicit 
time-dependency and works out as a Lagrange multiplier such as the total 
energy is an invariant of motion. Being the forcing mechanism not-positive 
defined, also the ideal viscosity will be not-positive defined. 
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The system now develops a net current and all the phenomenology of a 
dissipative physics. The natural question which arises is how much the orig- 
inal dynamics is preserved by this very-strong perturbation of the equation 
of motion. 

In JL7| some numerical simulation of eqs. ©-© have been performed 
showing that some of the main signature of the original physics are still 
present in the modified model with the advantage, in the latter, that also 
some analytical investigation can be carried out. In particular the most im- 
portant consequence of the chaotic hypothesis is the fluctuation theorem that 
is an exact parameterless prediction. This theorem concern the probability 
distribution function of the contraction rate of of the attractor surface ele- 
ment (for a detailed discussion see ||). Let us note that is very difficult, 
save for very particular systems, to test this prediction as the attracting sets 
usually are unknown. 

In the next section we will investigate a similar problematic in a class of 
dynamical systems for turbulent flows, called Shell- Models. In particular, we 
want to understand how much freedom is allowed in the choice of a viscous- 
modelization without perturbing too much the main physical framework, 
and/or quantifying the aspects of the perturbation, eventually. 
The goal consists in having a reversible dynamics showing in some limit (to 
be defined) the same physics of a turbulent dissipative flow. 

3 Time-reversible Shell Models 

One of the most intriguing problems in 3 dimensional turbulence is related 
to the understanding of the non linear dynamical mechanism triggering and 
supporting the energy cascade from large to small scales. Following the 
Richardson scenario that energy should be transferred downwards in scales, 
Kolmogorov Jl8j (K41) postulated that the energy cascade should follow a 
self-similar and homogeneous process entirely dependent on the energy trans- 
fer rate, e. This idea, with the assumption of local isotropy and universality 
of the small scales, eventually led to a precise prediction: 

S p (l) ee ((6v(l)Y) = C p ((e(l)y/3)P/* (5) 

where C p are constants and e(l) is the coarse-grained energy dissipation over 
a scale / supposed to be in the . inertial range, i.e. much smaller than 
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the integral scale and much larger than the viscous dissipation cutoff. If 
S p (l) ~ l<® and (e p (l)) ~ then 



C(p) = P/3 + r(p/3) (6) 

In K41 the e(Z) statistic is assumed to be /-independent, or r{p) = 0, 
implying ((p) = p/3, Vp, in particular £(2) = 2/3 or, equivalently, the energy 
spectrum going as A; _5//3 . On the other hand, there are many experimen- 



tal and numerical |19| , |20| , p2| results telling us that K41 scenario for 
homogeneous and isotropic turbulence is quantitatively wronged. Strong in- 
termittent bursts in the energy transfer have been observed and non trivial 
r(p) set of exponents measured. 

Shell models have demonstrated to be very useful for the understanding of 
many properties connected to the non-linear turbulent energy transfer |23| - 
[ |30|j . The most popular shell model, the Gledzer-Ohkitani-Yamada (GOY) 



model ( ||23|| -[p0[|), has been shown to predict scaling properties for ((p) (for 
a suitable choice of the parameters) similar to what is found experimentally. 

The GOY model can be seen as a severe truncation of the Navier-Stokes 
equations: it retains only one complex mode representative of all 

Fourier modes in the shell of wave numbers k between k n = k§2 n and k n+ \. 

Dynamical equations have the same qualitative structure of Navier-Stokes 
eqs., namely: 

— u n = N n [u) - vklu n + f n (7) 

where N [u] are the inertial nonlinear terms (see below) , while v is the molec- 
ular viscosity and f n a suitable forcing term acting only at large scale intro- 
duced in order to reach a (statistical) stationary state. 

The choice of the nonlinear term is dictated from the "locality assump- 
tion", i.e. only couplings with the nearest and next nearest shells are kept. 
In details the final eqs. are: 

— U n = ik n (u n+1 U n+2 + bu n+l U n _ 1 + Ctt n _ 1 tt n _ 2 ) — v ^a u n + ^n,n f, (8) 

where the the external forcing acts on a large scale n and b, c are two 
free parameters, but with the constraint 1 + 2b + 4c = 0, used for changing 
the "dimensionality" of the system |28|, |3(J, a popular choice which leads to 
results close to the 3d turbulent phenomenology is b = — 1/4, c = —1/8. Let 
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us stress that this choice of the parameters corresponds to have both energy 
and helicity conservation for a shell scale ratio equal to 2 and that whenever 
one has these two invariants, for any shell ratio, one has anomalous scaling 
exponents. 

At fixed molecular viscosity, v, the model develops a chaotic energy transfer 
to the small scales, with intermittent burst and deviation from K41 in good 
qualitative and quantitative agreement with what observed in true turbulent 
flow. 

The natural question which we would like to analyze in this paper is 
whether a reversible system obtained from the original dissipative equation 
by imposing a global constraint will allow us to reproduce the standard results 
and whether one can learn something of more about the strong-dissipative 
and far-from equilibrium structure of the stationary statistics. 

Therefore, following Gallavotti's suggestion we introduce a Lagrange mul- 
tiplier a[u] such as the eqs. of motion ([8]) preserve the total energy for any 
time, namely: 

— u n = N n [u] - a[u\k 2 n u n + f n . (9) 
In order to have the total energy E = 1/2 X] \u n \ 2 constant, one has to impose: 

«[ u ] = F k 2\ v 12 ( 10 ) 

Let us comment that in || |J different versions of reversible - hydrody- 
namical equations for a flow have been proposed depending on which macro- 
scopic observables one fixes by mean of the lagrangian multiplier. For ex- 
ample, equations similar to (|9]) but with a[u] chosen such that the total 
energy dissipation is conserved could in principle be used as well. In our 
view, guided from the phenomenological behaviour of turbulent flows, we 
believe that the only realistic constraint one can safely impose to the equa- 
tion is on the total energy. Constraining the total energy dissipation would 
put too much weight on the small scales statistics and would kill one of the 
most remarkable signature of turbulent flows: multifractal nature of energy 
dissipation. 

The goal of our study is to understand how the system move away from 
the stationary and equilibrium state that one obtain when f n = as soon as 
some energy pump and (perfect) energy sink are switched on (/„ > 0). 
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4 Numerical Results 



We first performed a benchmark numerical integration of a standard irre- 
versible and dissipative GOY model with fixed viscosity and forcing. This 
integration allows us to fix "physical realistic" values for the observables of 
the reversible dynamics. Numerical evolution was given by a forth-order 
Runge-Kutta algorithm, for a GOY model with N = 23 shells and a con- 
stant forcing on the first shell, / = 5 x (1 The integration time 
was several hundreds characteristic eddy turn-over times. We measured the 
structure functions and the average energy of the system. Afterwards, we 
integrated the reversible dynamical system keeping the total energy fixed to 
the mean value of the benchmark run. We kept all the other parameters of 
the model equal to those of the benchmark run, except for the value of the 
forcing which we let vary in order to switch continuously from a conservative 
equilibrium dynamics (/ = 0) to a conservative non-equilibrium dynamics 
(/>0). 

In fig. 1 we show the behaviour of functions: 



y 



U n U n+1 U n+2 + -U n -lU n U n+2 



P/3\ 

), (11) 



for p = 2 and for different values of the external forcing. E njP represent the 
power p/3 of energy flux from from shell nth to shell (n + l)th divided by 
k n . It is easy to show that in the inertial range the energy flux must be 



constant [27]. Therefore, it is natural to quantify the statistical properties of 
the energy transfer by measuring the scaling properties: S„ jP ~ k~^ p . One 
can observe that in the limit with vanishing forcing there is equipartition 
between degrees of freedom. In this case the viscous term is very low - the 
viscosity is proportional to the value of the forcing - and we have essentially 
a truncated-Euler system with gaussian probability distributions of the shell 
variables u n . When the forcing is increased there appears two different scal- 
ing ranges. In the first range (small fc's) it is clearly distinguishable an energy 
cascade. In the second range (large k's) the energy is in equipartition among 
the degrees of freedom. Likewise the probability distribution functions of 
shells in equipartition have all the same functional non-gaussian form. 
The range in which the energy cascade is observed is longer for higher forc- 
ing up to a critical forcing where the cascade range coincide with the inertial 
range of the original GOY model. 
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For even higher forcing the system falls in a stable fixed point in which all 
the energy is concentrated in the first shell and all the other shell variables 
are zero. 

We have checked that the cascade range is not due to finite size effects by 
performing a simulation with a larger number of shells (N = 28) and keeping 
all the others parameters constant. 

The same behaviour has been obtained in a model in which the dissipation 
term has been put only in the dissipative range (the last 7 shells). 

Likewise, we measured the scaling exponents of £ njP in the energy-cascade 
range, using the extended self similarity (ESS) in order to have more accu- 
rate fits. In the ESS one measures the behaviour of the structure functions 
of order p versus the structure function of order 3. In this way the scaling 
range is longer and measures of scaling exponents are more precise [ED] . 
In fig. 2 we plotted the function of order 6 and in fig. 3 we show the be- 
haviour of the scaling exponents compared to that of the GOY model. There 
is a clear intermittent energy cascade. 

On heuristic grounds, Gallavotti made the conjecture that a dissipative 
system and its time-reversible analog should be equivalent (in a statistical 
sense) if the time scale by which the viscosity reaches its mean value is shorter 
than the hydrodynamic time scales (i.e. the eddy turnover times). In this 
case the viscosity time-evolution would be confused with its average. 
Although we did not find a quantitative agreement with the conjecture of 
Gallavotti the qualitative behaviour is in the right direction. We defined the 
characteristic time r of a[u] as the decaying time of the autocorrelation of 
its logarithm, that is defined by: 



where x = log(a). We considered the logarithm of the viscosity function 
because the function itself has very large fluctuations (several order of mag- 
nitude larger than its average) and consequently the time average of its square 
has very long convergence times. We have found that in the case with the 
longer cascade range, i. e. in the system with statistical properties closer to 
the original GOY model, this characteristic time is shorter (see fig. 4) and 
consequently the time by which the viscosity reaches its average is shorter. 




(12) 
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Moreover let us stress that in this case a[u] has a smaller mean square value, 
i.e. smaller fluctuations. 

5 Conclusions 

The statistical mechanics of fully developed turbulence has features which 
are very different from those of the usual equilibrium statistical mechanics of 
Hamiltonian systems. Recently, Gallavotti proposed a dynamical ensemble 
equivalence between hydrodynamical dissipative system, e.g. the Navier- 
Stokes equations, and time-reversible systems. 

In this paper, we introduce and study a time-reversible dynamical system 
obtained, from a shell model for turbulence, changing the viscous term such 
as the energy is conserved. At small forcing values the system has statistical 
behaviours very close to those of a gas, i.e. energy equipartition and gaussian 
statistics. At increasing the forcing one has a non-equilibrium statistical 
stationary state with an energy cascade and anomalous scaling laws similar 
to those observed in turbulence. 

The dynamical ensemble equivalence seems to be satisfied at least for the 
typical observables measurable in turbulent flows. 

A relevant open problem remains for defining precisely the class of con- 
straints allowed for the Dynamical Ensemble Equivalence hypothesis to hold. 

6 Acknowledgement 

This work was supported by INFN (Iniziativa Specifica Meccanica Statistica 
Fill) and INFM (Progetto di Ricerca Avanzato TURBO). 
We thank F. Bonetto and G. Gallavotti for many stimulating discussions and 
continuous encouragement. We also thank L.P. Kadanoff and R. Pasmanter 
for useful discussions. 



13 



References 

[1] U. Frisch, Turbulence, (Cambridge University Press, Cambridge, UK, 
1995). 

[2] J. Sommeria, J. Fluid Mech., 170, 139 (1986). 

[3] C. Beck and F. Schlogl, Thermodynamics of chaotic systems, (Cam- 
bridge University Press, Cambridge, UK, 1993). 

[4] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995). 

[5] G. Gallavotti , Phys. Lett. A 223, 91 (1996). 

[6] G. Gallavotti, Dynamical Ensembles Equivalence in Fluid Mechanics, 
preprint in chao-dyn@xyz.lanl.gov, (j 9605006. 

[7] Z. S. She and E. Jackson, Phys. Rev. Lett 70, 1255 (1993). 

[8] R. H. Kraichnan, J. Atmos. Sci. 33, 1521 (1976). 

[9] M. Briscolini and P. Santangelo, J. Fluid Mech. 270, 199 (1994). 

[10] T. Bohr, M. H. Jensen, G. Paladin and A. Vulpiani, Dynamical Systems 
Approach to Turbulence, (Cambridge University Press, Cambridge, UK, 
1997). 

[11] L. Kadanoff, Physics Today, September 11 (1995). 

[12] R. H. Kraichnan and D. Montogomery, Reports Prog. Phys. 43, 547 
(1980). 

[13] R. H. Kraichnan, J. Fluid Mech. 59, 745 (1973). 

[14] J. Sommeria, C. Nore, T. Dumont and R. Robert, C. R. Acad. Sci. 
Pans, serie 7/312, 999 (1991). 

[15] J. Miller, Phys. Rev. Lett. 65, 2137 (1990). 

[16] R. A. Pasmanter, Phys. Fluid, 6, 1236 (1994). 



14 



[17] F. Bonetto, G. Gallavotti and P. Garrido, Chaotic Principle: an exper- 
imental test preprint in mp_arc@math.utexas.edu, ft 96-182. 

[18] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 9 (1941). 

[19] C. Meneveau and K. R. Sreenivasan, J. Fluid Mech. 224, 429 (1991). 

[20] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, C. Massaioli & S. 
Succi, Phys. Rev. E 48, R29 (1993). 

[21] R. Benzi, S. Ciliberto, C. Baudet and G. R. Chavarria, Physica D, 80, 
385 (1994). 

[22] R. M. Kerr, J. Fluid Mech., 211, 309 (1990). 
[23] E.B. Gledzer, Sov. Phys. Dokl. 18, 216 (1973). 

[24] M. Yamada and K. Ohkitani, Prog. Theor. Phys. 81, 329 (1989); J. 
Phys. Soc. Jpn. 56, 4210 (1987); Phys. Rev. Lett 60, 983 (1988); Prog. 
Theor. Phys. 79, 1265 (1988). 

[25] M.H. Jensen, G. Paladin and A. Vulpiani; Phys. Rev. A 43, 798 (1991). 

[26] L. Biferale, A. Lambert, R. Lima and G. Paladin; Physica D 80 105 

(1995) . 

[27] D. Pisarenko, L. Biferale, D. Courvoisier, U. Frisch and M. Vergassola, 
Phys. Fluids A65 2533 (1993). 

[28] L. Kadanoff, D. Lohse, J. Wang and R. Benzi; Phys. Fluids 7 617 (1995). 

[29] R. Benzi, L. Biferale, R. Kerr and E. Trovatore, Phys. Rev. E 53 3541 

(1996) . 

[30] L. Biferale and R. Kerr, Phys. Rev. E 52, 6113 (1995). 



15 



FIGURE CAPTIONS 



FIGURE 1: 

S n ,2 vs k n in log-log scale obtained from an integration of the model 
with N = 23, Kq = 6.25 x 10~ 2 and different values of the forcing: / = 
5 x 1(T 3 (1 +i) (plus), / = 6 x 1(T 3 (1 +i) (squares), / = 8 x 10" 3 (1 +i) 
(cross). Diamonds represent the results obtained in the benchmark 
integration (i.e., the original GOY model). 

FIGURE 2: 

S n 6 vs S„ 3 in log-log scale obtained for the system with the same 
parameters as in figure 1 and / = 8 x 10 _3 (1 +i) (diamonds). The K41 
line (dashed line) is also shown for comparison. 

FIGURE 3: 

( p vs p for the GOY model(diamonds) and for its reversible analog. 
Squares are obtained for / = 8 x 10~ 3 (1 + i) and pluses for / = 5 x 
1CT 3 (1 + i). The anomalous exponents have been calculated using the 
ESS. 

FIGURE 4: 

C(t) vs t for two different values of the forcing: / = 5 x 10~ 3 (1 + i) 
(solid line) and / = 8 x 1CT 3 (1 + i) (dashed line). 
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